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ABSTRACT 

We present tentative evidence for primordial clustering, manifested as an excess of color-selected 
objects in the field of the QSO SDSS J 1030+0524 at redshift z = 6.28. We have selected objects red 
in «775 — Z850 on the basis of Hubble Space Telescope Advanced Camera for Surveys imaging of a field 
centered on the QSO. Compared to data at comparable depth obtained by the GOODS survey, we 
find an excess of objects with (1775 — zsso) > 1.5 in the QSO field. The significance of the detection 
is estimated to be ^97% on the basis of the counts alone and increases to 99.4% if one takes into 
account the color distribution. If confirmed this would represent the highest redshift example of galaxy 
clustering and would have implications on models for the growth of structure. Bias-driven clustering 
of first luminous objects forming in the highest peaks of the primordial density field is expected in 
most models of early structure formation. The redshift of one of the candidates has been found to be 
z = 5.970 by our spectroscopy with Keck I/LRIS, confirming the validity of our color selection. 
Subject headings: cosmology: early universe; galaxies: formation; high-redshift 
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crowave background (Kogut et al. 2003) would place the 
(initial?) reionization at z ~ 10-20. 

A simple picture of a clean-cut, universally synchro- 
nized Gunn-Peterson trough now appears unlikely. The 
key issue is the dumpiness of the IGM, and the gradual 
development and clumpy distribution of the first ioniz- 
ing sources, either protogalaxies or early AGN (see, e.g., 
Miralda-Escude et al. 2000). Reionization is expected 
to occur gradually as the UV emissivity increases (cf. 
McDonald & Miralda-Escude 2001) with the lowest den- 
sity regions becoming fully reionized first. This is also 
suggested by modern numerical simulations (e.g., Gnedin 
& Ostriker 1997; Ciardi et al. 2003; Gnedin 2004) which 
predict an extended period of reionization, starting at 
z ~ 15 or even higher and ending at z ~ 6 or so (see also 
Wyithe & Loeb 2003, Cen 2003, Haiman & Holder 2003, 
Somerville et al. 2003). 

Large samples of normal galaxies are now studied out 
to z ~ 4.5 (Steidel et al. 1999, Dawson et al. 2004). Sur- 
veys like, e.g., GOODS (Giavalisco et al. 2004, Dickinson 
et al. 2004) and HUDF (Beckwith et al. in preparation; 
see also Bunker et al. 2004, Yan & Windhorst 2004), have 
provided us with the first significant samples of galaxies 
at z > 5. However, if our understanding of galaxy forma- 
tion thrmich hierarchical chisterinp' is rnrrcrt wp pynprt 



1. INTRODUCTION 

Exploration of the early structure formation and the 
reionization era is now one of the most active areas of 
cosmology; see, e.g., Madau (2000), Loeb & Barkana 
(2001), Barkana & Loeb (2001), Stiavelli (2002), or Djor- 
govski (2004) for recent reviews and references. The first 
hints of an approach to reionization, signaled by the dra- 
matic increase in the intergalactic hydrogen absorption 
at z ~ 6, was detected in the spectra of high-z QSOs 
(Becker et al. 2001, Djorgovski et al. 2001, White et al. 
2003, etc.). In particular, there was the possible detec- 
tion of an extended, optically thick absorption due to 
neutral hydrogen at \ re st < 1216 A (Gunn & Peter- 
son 1965) in the spectrum of QSO SDSS J1030+0524 at 
z — 6.28 (Fan et al. 2001). These results strongly suggest 
the end of the cosmic reionization at z w 6. However, 
WMAP's tentative detection of a large amplitude signal 
in the temperature-polarization maps of the cosmic mi- 
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tcring at redshift around 6 has been found in the UDF 
(Malhotra et al. 2004, see also Ouchi et al. 2005). This 
leads us naturally to the study of z > 6 QSO environ- 
ments as they should be the most clustered environments 
at these very high redshifts. 

Indeed, a generic expectation in most models of galaxy 
formation is that the most massive density peaks in the 
early universe are likely to be strongly clustered (Kaiser 
1984; Efstathiou & Rees 1988). The evidence for such 
bias is already seen in large samples of Lyman-break 
galaxies at z ~ 3-3.5 (Steidel et al. 2003), and in Lyman 
a selected galaxy samples (e.g., Venemans et al. 2003); it 



2 



Stiavelli ct al. 



should be even stronger at higher redshifts. An "excess" 
in the number of galaxies and in the density of star for- 
mation was also discovered in a systematic Keck survey 
of fields centered on known z > 4 quasars (e.g., Djorgov- 
ski 1999, Djorgovski et al. 1999, 2003, and in prep.). The 
high metallicity associated with QSOs — even at z > 6 — 
is often interpreted as evidence that they are located at 
the center of massive (proto-)galaxies, thereby corrobo- 
rating the overall picture. These arguments justify the 
expectation that QSOs at z ~ 6 most likely highlight 
some of the first perturbations that become non-linear 
in the density distribution of matter. 

This paper is devoted to the study of the environment 
of the QSO SDSS J1030+0524. Section 2 describes the 
data and our data reduction techniques. Section 3 de- 
scribes the catalogs and the candidate list. Section 4 
presents a spectrum obtained for one of our candidates. 
Section 5 presents and discusses our results. 

2. DATA REDUCTION AND ANALYSIS 

We have observed a field near SDSS J1030+0524 us- 
ing the Advanced Camera for Surveys (ACS) onboard 
the Hubble Space Telescope. The Wide Field Camera 
(WFC) was used to secure an integration of 5,840 s in 
F775W (hereafter i 775 ) and of 11,330 s in F850LP (here- 
after Z85o). The relative exposure times in these filters 
have been chosen so as to have a depth in these two filters 
comparable to that of the GOODS survey (Giavalisco 
et al. 2004) allowing us to use GOODS as a reference 
field sample. Since SDSS J1030+0524 is about 30 arcsec 
from a 12th magnitude star, we placed the star in the in- 
terchip gap of the WFC. This avoids any bleeding along 
the CCD columns. The images were dithered to improve 
the PSF sampling and we obtained a total of 7 images in 
F775W and 9 images in F850LP. 

The data were processed by the ACS pipeline 
(CALACS), which performs the basic data reduction 
steps of bias subtraction, dark subtraction and flat field- 
ing of the raw data, resulting in the flat-fielded images 
(the fit files) . The post-pipeline processing steps include 
the creation of individual weight maps for each fit file, 
drizzling to the same scale as GOODS (0.03 arcsec/pixel) 
using MultiDrizzle (Koekemoer et al. 2002), correct- 
ing for correlated noise introduced by MultiDrizzle, and 
masking any bad pixels and pixels that were less than 
10% of the total exposure time. The final MultiDrizzle 
run was performed with parameters pixscale = 0.6 and 
pixfrac = 1.0 and with final _wht_type = ivm to make use 
of individual weight maps. Special care was taken in pro- 
ducing the weight maps used as an input to MultiDrizzle, 
which also helps eliminate many spurious detections from 
the resulting object catalogs. 

Source detection was done on the zsso-band images, 
using the same version of Sextractor (Bertin & Arnout 
1996) used by the GOODS collaboration and the same 
input parameters in order to facilitate a straightforward 
comparison between the two catalogs. The software was 
run with the MultiDrizzle output image (drz file) and 
weight map (converted to an rms map). The Sextrac- 
tor parameters used attempt to maximize the number 
of faint sources found while minimizing the number of 
false detections. We performed isophotal photometry, 
mag_auto photometry, photometry with a 16 pixel fixed 
diameter and photometry at the half light radius. The 



photometric errors calculated by Sextractor are mea- 
sured from the weight maps and take into account the 
background and instrumental noise and exposure time 
information. We verified that the various photometric 
measurement were consistent and for the rest of the anal- 
ysis we used the mag_auto photometry. We have only 
considered sources detected in Zsso a t S/N > 5. For these 
sources we have computed the (4775 — zsso) color directly 
from the 1775 magnitude when the 1775 flux was measured 
to better than 2a or otherwise from the 2a upper limit to 
the 1775 flux. We have run completeness simulations in 
the SDSS J1030+0524 field and found that our photom- 
etry has comparable completeness to GOODS. We find 
that we reach 50% completeness for compact galaxies 
at z 8 5o — 27. Finally, we have verified that our cata- 
log contains no contamination by noise peaks mistakenly 
identified as galaxies by running a negative image test 
(Dickinson et al. 2004, Yan & Windhorst 2004). 

3. CANDIDATE OBJECTS ATz = 6 

The final Sextractor catalog contains 1551 objects. 
Fig. la shows the counts for the SDSS J1030+0524 field 
and those for the GOODS N and S fields corrected to the 
same area. The QSO is recovered as the reddest object 
in the field. Fig. lb shows the color distribution of all 
objects with z < 26.5. The basic distribution in both 
counts and color are similar in the SDSS J1030+0524 
area and in GOODS. This is an independent test that 
the GOODS catalog and ours are generally consistent. 

Excluding the QSO, a total of 8 objects have (1775 — 
Z85o) > 1-3, S/N > 5 and ZS50 < 26.5. We expect that a 
high fraction of galaxies with (1775 — zsso) > 1-3 will be 
at z > 5.5 as Malhotra et al. (2005) have shown that at 
^850 < 27 about 90% of these objects are truly at high- 
rcdshift. At our color cut, and in agreement with Willott 
et al. (2004), we do not find a significant overall source 
density excess between the SDSS J1030+0524 field and 
GOODS. However, when we adopt the more conservative 
cut (1775 — Z85o) > 1-5, wc still find 7 objects (excluding 
the QSO) in the SDSS J1030+0524 field against the 3.33 
expected from GOODS. This is illustrated graphically in 
Fig. lc where we show the counts of objects redder than 
(*775 — ^85o) = 1-5. The objects are generally resolved as 
shown in Fig. 2 for two objects from this selection. 

4. A FIRST ATTEMPT AT SPECTROSCOPIC 
CONFIRMATION 

Spectroscopic observations of eight of the color-selected 
candidates were obtained using the LRIS instrument 
(Oke et al. 1995) on the W.M. Keck Observatory's 10-m 
Keck-I telescope in multislit mode. The data were ob- 
tained on the nights of UT 2004 December 10 and 15, 
in moderate to good conditions, with an effective total 
exposure time of 2.8 hours per night. All observations 
employed a 400 lines mm -1 grating (Ablaze = 8500 A), 
giving an effective resolution of FWHM ~ 8 A. The 
data were reduced using standard procedures and flux- 
calibrated using archival data obtained in photometric 
conditions through the same instrument configuration. 

Only one of the targeted sources, J103024.08+ 
052420.4, shows an emission line, allowing an unambigu- 
ous redshift identification. Its spectrum, derived from 
the UT 2004 December 10 data when conditions were 
more favorable, is shown in Fig. 3. The emission line 
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reproduced on UT 2004 December 15. The spectrum 
shows a single, isolated, asymmetric, high equivalent 
width emission line, characteristic of high-redshift Lya 
emission (cf. Stern et al. 2005). A weak continuum decre- 
ment across the emission line is also evident, character- 
istic of absorption due to the Lya forest at high redshift. 
We derive a redshift of z = 5.970 for this source, con- 
firming that our i-drop criteria reliably identify z ~ 6 
galaxies. 

The spectra of the remaining targets do not have a 
sufficient signal-to-noise to make any other firm identifi- 
cations at this point. Further analysis of these and other, 
related data will be presented elsewhere. 

5. DISCUSSION AND CONCLUSIONS 

The probability of finding 7 objects when one expects 
3.33 is 0.7% if the distribution is Poissonian. However, 
the true probability would probably be higher because of 
clustering. Indeed, by randomly overlaying WFC fields 
within GOODS we find that the probability of one such 
field containing 7 objects with (1775 — zsso) > 1-5 is about 
3%. 

Luckily we do not only have the number of red objects 
but also their distribution in color and this also turns 
out to be different from that within GOODS. In Fig. Id 
we show the cumulative distribution of objects within 
the SDSS J1030+0524 field and in GOODS (renormal- 
izcd to the same area). Clearly, the SDSS J1030+0524 
object distribution appears to be different from that of 
GOODS. We have applied a Kolmogorov-Smirnov test 
to determine whether the distribution of the 8 sources in 
SDSS J1030+0524 redder than {i 775 - z 850 ) = 1.3 and 
that of the 252 in GOODS with the same color cut are 
drawn from the same parent distribution. The result is 
a difference value of D = 1.72 corresponding to a proba- 
bility 0.55%. 

Thus, we find that the SDSS J1030+0524 field appears 
to have an excess of red objects. Our result is in agree- 
ment with Willott et al. (2004) as our data set is deeper 
and the majority of the excess sources are fainter than 
^850 = 25.5. 

If the SDSS J1030+0524 does indeed have an excess 
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of red sources, one may wonder whether their redshift 
is the same as the QSO. The selection (1775 — zsso) > 
1.5 should be effective at selecting objects at redshift 
greater than 5.9. Indeed, the one object for which we 
could obtain a spectroscopic confirmation satisfies this 
constraint. Most of the sources are too faint to be able 
to measure a color as red as that of the QSO but 3 of the 7 
sources are detected as dropouts, i.e. their 1775 band flux 
is measured at less than 2a. It should also be noted that 
SDSS QSOs at z > 6 are extremely rare objects and are 
presumably associated to rare density peaks which may 
correspond to perturbations on very large scales. The 
size of large scale structures seen at lower redshift can go 
up to at least 50 comoving Mpc which would correspond 
to Sz — 0.12 at z = 6.28. Thus, we expect some, but 
not all, of the excess sources to be at z — 6.28. We do 
not know whether J103024.08+052420.4 belongs to the 
same overdensity as the QSO as the redshift difference 
is sizeable. However, one can speculate that if it did, it 
would imply a minimum size for the overdensity of ~ 130 
comoving Mpc along the line of sight. If the overdensity 
had this typical size also in the transverse direction, this 
would translate into an angular extension of ~ 50 arcmin. 

The discovery of a density excess near SDSS 
J1030+0524 is potentially very interesting and requires a 
further confirmation. Additional spectroscopic follow-up 
is now in progress. The presence of a very extended over- 
density could be verified by repeating this search over a 
much larger area. We also plan to repeat this study on 
4 more fields around redshift ~ 6 QSOs, which will be 
reported elsewhere. 
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Fig. 1.— Comparison of sources in the field around SDSS J1030+0524 (thick lines) with those in GOODS (thin lines). The GOODS 
counts have been renormalized by the ratios of the area of the two surveys and only sources with S/N > 5 are shown. Panel (a) shows 
the total counts. All remaining panels include only objects brighter than 2850 = 26.5. Panel (b) shows the color distribution within the 
two catalogs. Both panels (a) and (b) illustrate how the generic properties of objects in the two catalogs are comparable. Panel (c) shows 
the counts of objects redder then 1775 — 2850 = 1-5. Panel (d) shows the number of objects redder than a given 1775 — 2:850 as a function 
of 1775 — z 850- I n panel (d) we have omitted 4 objects that didn't look convincing by visual inspection (no similar check was done on the 
GOODS sources). The field of SDSS J1030+0524 appears to contain an excess of sources with 1775 — 2350 > 1.5. 
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Fig. 2.— Images of two of the reddest objects in the SDSS J1030+0524 field: J103021. 74+052510.8 (top) and J103024.08+052420.4 
(bottom). For each we show the image in F775W on the left and the one in F850LP on the right. Both objects appear resolved in the 
F850LP images. 
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Fig. 3.— Spectrum of J103024. 08+052420. 4 at z = 5.970, obtained with LRIS on the Keck-I telescope. The total exposure time is 2.8 hr, 
and the spectrum was extracted using a 1'.'2 X l'/5 aperture. The spectrum has been smoothed using a 5.5 A boxcar filter. The dotted 
spectrum in the main panel indicates the corresponding la error spectrum, assuming Poisson fluctuations of the sky spectrum. 



